pedmf boundary conditions #4106
Merged
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Purpose
This PR updates the treatment of boundary conditions for prognostic EDMF:
1- Apply surface fluxes to both grid-mean and updraft fields equally.
2- Seed a small updraft area fraction in the first cell when the surface buoyancy flux is positive, allowing it to grow even if initially zero.
3- Instead of overwriting the first-cell area fraction, increase the entrainment rate when buoyancy flux is positive, ensuring the area reaches the prescribed value in one timestep. This ensures that tracer masses are conserved.
4- Instead of overwriting mse and q_tot, assume a near-surface entrainment in the first cell, with tracer values of the entraining fluid provided by the function
sgs_scalar_first_interior_bc(mean + perturbation).Content